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ABSTRACT 

The singular bifurcations in a dispersive billiard are discussed in terms of symbolic 
dynamics and is compared to an example of a bifurcation tree in a smooth potential. 
Possible generalizations to other smooth potentials are discussed. 
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1 Introduction 



The aim of this article is to describe the bifurcation structure of some families of 2- 
dimensional Hamiltonian systems. Bifurcations in one dimensional unimodal maps 

are well understood, and the Metropolis, Stein, Stein (MSS)0 theory uses symbolic 
dynamics to establish that the order of bifurcations is universal for all unimodal 
maps. However, there is no general theory of bifurcations in Hamiltonian systems 
analogous to the MSS theory for unimodal maps. The method presented here 
predicts the bifurcations of a periodic orbit, including all periodic orbits emanating 
from a bifurcation, and the non-periodic (chaotic) orbits born at the accumulation 
points of bifurcations. We refer to these orbits as the members of one bifurcation 
family, and we show that the members of a family are well defined and do not change 
with small changes of the Hamiltonian. However, the ordering in the parameter 
space of the bifurcations within a family may be different for different smooth 
potentials, and the ordering of different families in the parameter space is not known 
a priori. 

We are not able to define rigorously the class of smooth Hamiltonians to which 
the method described here applies, but the Hamiltonian must have some similarities 
with a billiard system. This includes physical interesting examples such as the 

classical collinear CO2 molecule0 ; and possibly the classical collinear He atom0 

and the St0rmer problem^. The smooth potential we use as an example is (for 

reference see E' 0) 



where variation of the parameter a gives rise to a family of Hamiltonians. In the 
limit a — > the system reduces to the hyperbola billiard and all orbits can be 

assigned a well defined symbolic dynamics^' ^. As a increases, orbits disappear 
through inverse bifurcations. The bifurcations are period doublings, symmetry 
breaking, tangent bifurcations and bifurcations of higher order, and they always 
imply stable orbits; hyperbolic orbits become elliptic orbits for some parameter 

interval and then disappear 0, as shown in the examples below. 

A billiard system, e.g. the 3-disk billiard discussed below, with a parameter r, 
also has a bifurcation scenario. For a parameter value r c , an infinite number of 
orbits are pruned, but in the billiards considered here, no orbits become stable. 

We begin by reviewing the bifurcations in one dimensional unimodal maps. We 
then describe the simplest example in a 3-disk billiard, and discuss in detail a 
bifurcation sequence in the 4-disk billiard, and compare it with the bifurcations 
found by Dahlqvist and Russberg for the Hamiltonian ([I]). 




(1) 
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2 Unimodal maps 



Here we motivate our approach by showing that bifurcations in a one-dimensional 
logistic map are analogous to bifurcations in a smooth potential, and that bifurca- 
tions in the tent map are analogous to bifurcations in a billiard. The logistic map 
is defined as 

x t+ i = ax t (l - x t ) 

and the symbolic dynamics we use is symbol 1 if x > 1/2 and symbol if x < 1/2. 
The MSS bifurcation scenario starts out as follows: at a = 1 the fixed point 1 
becomes stable, at a = 3 the fixed point 1 becomes unstable while the period 2 
orbit 10 is born as a stable orbit, at a = 3.4494 . . . the orbit 10 becomes unstable 
and the period 4 orbit 1011 is born as a stable orbit, and so on. The MSS theory 
predicts period doubling of an orbit into a periodic orbit of twice the length. The 
symbolic dynamics description of the new orbit is obtained by writing the original 
symbol sequence twice and then changing the last symbol from to 1 or from 1 
to 0; orbit 1011 period doubles into 10111010, etc. At the accumulation point 
a = 3.5714 ... all such orbits exist and are unstable. We call this collection of 
orbits the period doubling family of the fixed point 1. In the same way one finds 
the period doubling family of the period 3 orbit 100, and other orbits born by 
tangent bifurcations. To be more precise, we can identify symbol with x < 1/2 
and symbol 1 with x > 1/2 for stable orbits in this family, only when the orbit 
has passed through the superstable point x — 1/2. Between the bifurcation point, 
where the orbit is born, and the superstable point, the naive symbolic description 
of the stable orbit is the same as for a different unstable orbit. This unstable orbit 
is either an orbit of the same length created in a tangent bifurcation, or the orbit 
which became unstable at a period doubling bifurcation. 
In the tent map 



this bifurcation scenario is slightly altered. The fixed point 1 exists only when 
r > 1. Since the slope at this fixed point is — r < — 1, the fixed point is unstable 
and the period 2 orbit 10 exists when r > 1. This orbit has stability — r 2 < — 1, 
and is also unstable. By induction all the orbits in the period doubling family of 1 
exist and are unstable for r > 1 and consequently r c = 1 is the critical parameter 
value where all orbits in this family bifurcate simultaneously. For r = r c all points 
in all orbits have x- values equal 1/2 and when r increases from r c the points in the 
orbits spread out on the x-axis. The critical value for the period 3 family 100 is 
r c = (1 + v / 5)/2 = 1.6180. . ., with all orbits in the period doubling family born 
unstable, and the same happens for all orbits born in tangent bifurcations. There 
may also be other orbits born at these critical parameter values. 




rx t x < 1/2 

r(l - x t ) x > 1/2 
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MSS universality is a consequence of the fact that the parameter space is one- 
dimensional for a unimodal map, and there is only one possible ordering of bifur- 
cations for any monotone variation of the parameter. An n-modal one- dimensional 
map has more complicated bifurcations in an n-dimensional parameter space. For 
higher dimensional maps, the situation is totally different; the parameter space is 
in principle always infinite dimensional. 



3 Bifurcations in disk billiards and potentials 

The billiards that we use here as our standard examples of dispersive billiards 
consist of a point particle scattering elastically between N circular disks in the 
plane. The basic symbolic description is obtained by enumerating the disks, and 
letting a bounce off disk number s at the integer time t be denoted by the symbol 
St- Other symbols can be made from St G {1, . . . , N}, the number of symbols can 

be reduced^' ^H, and a well ordered alphabet can be constructed S —I, but for 
simplicity we shall use only the basic symbols here. 

The simplest chaotic billiard consists of 3-disks of radius 1, with the distances 
between the centers of the disks given by parameters ri2, ri3 and r23- If the disks 
are sufficiently separatedE3 ; all orbits that can be described by any infinite symbol 



sequence . . . s t _is t s t+ i . . . exist, except orbits with two successive bounces off the 



same disk. As shown in ref. ^ E3, a bifurcation in the billiard takes place when 
a free path between two bounces is tangent to a disk, or when the bounce of the 
particle has an incoming and outgoing direction parallel to the edge of the disk. 
If the parameter r changes further the orbit is not admissible in the system; we 
call such orbit pruned. It is clear that at the critical bifurcation parameter value 
r c the orbit that is tangent to the disk is indistinguishable in configuration space 
from the orbit that bounces parallel to the edge of the disk. From this fact we can 
establish the set of orbits that bifurcate at r c and at the same point in phase space. 
Analogously to the one-dimensional tent map we call this set of orbits a bifurcation 
family. 

Consider as an example the periodic orbit 123132 drawn in figure [3] a) for r\i = 
ri3 = = 2.5. To obtain simple figures we have chosen the distance between the 
centers of the disks 2 and 3 to be the parameter r = r23 while keeping ri2 = ri3 = 
2.5. This orbit bifurcates at r c = 4.1044 . . ., as shown in figure [l] b) where the line 
from disk 2 to disk 3 is tangent to disk 1, and the orbit is not admissible for r > r c . 
Figure [I] c) shows the periodic orbit 12131312 where the particle bounces off disk 
1 instead of moving directly between disk 2 and disk 3. This orbit also bifurcates 
at r c = 4.1044. . . where, in the configuration space, it coincides with the orbit of 
figure |l] a) . 

Clearly all orbits that coincide with the orbit of figure [I] b) for r c bifurcate exactly 
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at this parameter value, and the symbolic description of these orbits is obtained 
by either including or not including the symbol for disk 1 each time the particle is 
tangent to disk 1. One example is when we chose to include the symbol 1 each time 
the particle goes from disk 3 to disk 2 but not when going from disk 2 to disk 3. 
This gives the periodic orbit 1231312 drawn in figure [I] d). This orbit has the same 
length in configuration space as the first orbits for the bifurcation parameter r c , but 
it has a broken symmetry and we can think of it as the orbit born at a symmetry 
breaking bifurcation in a smooth potential. Examples of other periodic orbits in the 
bifurcation family are 1231321231312, which does not bounce at the tangency point 
the first 3 times, and then bounces once, and the orbit 123131212131312 that does 
not bounce the first time, but then bounces the next 3 times at the tangency point. 
The length of the symbol sequence of these two orbits is approximately twice the 
length of the symbol sequence of the short orbits, and for the bifurcation parameter 
value r c , the length of these orbits in configuration space is exactly twice the length 
of the short orbits. These orbits could be born at period doubling bifurcations 
in a smooth potential. Orbits 3,4,5,. . .times the length of the original orbits are 
constructed the same way. 

All orbits in this bifurcation family are of the form 

...212s_i313s 212si313s 2 ... 

where Sj stands for a symbol 1 or no symbol. Periodic sequences of Sj correspond 
to periodic orbits while non-periodic sequences of Sj correspond to chaotic orbits. 

We observe that the scenario of bifurcations is similar to the tent map; all 
bifurcations take place at one parameter value r c . The number of orbits in a family 
here is however larger than in the one dimensional unimodal map period doubling 
scenario; in one dimension we repeated the symbol sequence twice and changed the 
last symbol, while here we repeat a short symbol sequence and then have different 
possibilities for changing a symbol. We expect the bifurcations in a smooth potential 
to be more complicated and with a many dimensional parameter space. 

To compare a billiard system with the smooth Hamiltonian ([J) we investigate the 
bifurcation family of the orbit 1(32) 4 4(23) 4 in the 4-disk billiard shown in figure^ a). 
The parameter in the 4-disk billiard is r = ri 2 = r 2 3 = = r^\. This orbit 
bifurcates when the line from disk 1 to disk 3 is tangent to disk 2, and because 
of the symmetry of the orbit, at the same parameter value the line from disk 4 
to disk 2 is tangent to disk 3. The orbit that bounces off each of the tangent 
points is the orbit 12(32) 4 343(23) 4 2 drawn in figure |2| b). We can now construct all 
possible combinations where the orbit does or does not bounce off the disk 2 or 3 at 
the tangent points. The shor test possibilities are: 1(32) 4 34(23) 4 2, 1(32) 4 343(23) 4 , 
1(32) 4 4(23) 4 2, 1(32) 4 43(23) 4 2 and 1(32) 4 343(23) 4 2, and orbits symmetric to these. 
The general form of this bifurcation family of orbits is 

. . . ls_i(32) 4 t_ 1 4t (23) 4 s lsi(32) 4 t 1 4t 2 (23) 4 s 2 ... (2) 
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with Si either 2 or no symbol, and t{ either 3 or no symbol. 

The position on the edge of disk 1 for the bounce, is shown in figure [3] for some 
of these orbits as a function of the parameter r. The critical parameter value is 
r c = 2.0312 . . . and the orbits are not physically admissible for r < r c . Figure [3] 
shows that when r > r c the positions spreads out with orbit 1(32) 4 4(23) 4 having 
the largest value and the orbit 12(32) 4 343(23) 4 2 having the smallest value. 

The bifurcation diagram for some of the orbits in this family for the Hamiltonian 

([!]) is described in ref.0 and is sketched in figure f|. Only the short orbits are 
included and the drawing is not correctly scaled in the parameter value, and the 
horizontal positions are chosen only for better visualization (for numerical values see 

figure 4 in ref.0). A dotted line is a hyperbolic orbit while a solid line is an elliptic 
orbit. The figure shows that the two orbits 1(32) 4 4(23) 4 and 12( 32) 4 343(23) 4 2 are 
born together, and that 1(32) 4 4(23) 4 is the hyperbolic orbit while 12(32) 4 343(23) 4 2 
is the elliptic orbit. The elliptic orbit then has bifurcations where it bifurcates 
into orbits belonging to the family (^j). This family structure of symbols for this 

Hamiltonian (JJ) is conjectured by Dahlqvist and RussbergB. By finding longer 
orbits we expect each elliptic branch of the tree to go through an accumulation of 
bifurcations similar to the logistic map. There is no reason to expect our choice of 
the parameter a to be in any way a universal choice, or that this should be a universal 
one-dimension parameter bifurcation scenario like for the unimodal map. A different 
parametrization of the Hamiltonian may give a slightly different bifurcation tree, 
but always including exactly the orbits in the family (0). The number of orbits in 
the bifurcation family is larger than the period doubling family in the logistic map, 
and this implies that some orbits have to bifurcate several times. For example, in 
the logistic map no orbits bifurcates more than once, while in the example above 
the orbit 12(32) 4 343(23) 4 2 bifurcates three times. 

The example shows that if we know the bifurcation family in a billiard system 
we can predict the expected bifurcations in a smooth potential and the geometrical 
shape of the orbits created in a bifurcation. We expect a potential with 3 or 4 
smooth hills to exhibit bifurcation trees of similar structure. The theory of bifur- 
cation in smooth potentialsGl' 113 predicts bifurcations when a eigenvalue A of 
the monodromy matrix has rational phase on the unit circle, or when it becomes 
real. The conservation of the Poincare index gives some restrictions on the possi- 
ble bifurcations. The theory does however not predict the geometrical structure of 
the bifurcating orbits, nor tells anything about how many times an orbit bifurcates 
with the same eigenvalue. Study of the billiard bifurcations adds this information 
to the theory. This may give a unique symbolic description of all unstable and 
stable orbits in a potential. 

An important implication of the billiard studies is that if the structure of bifur- 
cations in the 4-disk billiard and in the smooth potential are similar, then we expect 
to be able to describe the allowed orbits in the smooth potential by a pruning front 
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similar to the pruning front for the billiard"' I— U. This fives a method to find an 
approximation of a Markov partition of the system which is useful in calculations 
of thermodynamical quantities 



4 Conclusions 

We have here illustrated the utility of comparing bifurcations in billiard systems 
and in smooth Hamiltonian systems by an example comparing disk billiards and 
the (x 2 y 2 ) l / a potential. The bifurcations is easy to find in billiards and we obtain 
the pruning front which determines the forbidden orbits and the symbolic dynamics 
for the system. It is more difficult to find the bifurcations in a smooth potential 
and the existence of a pruning front in that case is not yet established. We have 
shown here that by using symbolic dynamics we can find bifurcation families of 
orbits and that these families have the same structure in the billiard and in the 
smooth potential. This indicates that the pruning front is similar for the billiards 
and for smooth potentials. Knowledge of the bifurcations in billiards can also be of 
help in numerical searches for bifurcations in a potential. 
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1. Examples of orbits that belongs to the same bifurcation family for the 3 disk 
billiard, a) The orbit 123132 with r 12 = r 13 = r 23 = 2.5. b) The orbit 123132 
with r 12 = r 13 = 2.5 and r 23 = r c = 4.1044 . . .; the bifurcation point for the 
orbit, c) The orbit 12131312, r 12 = r 13 = r 23 = 2.5. d) The orbit 1231312, 
rw = r 13 = r 23 = 2.5. 

2. Examples of orbits that belong to the same bifurcation family in the 4 disk bil- 
liard with r = 2.5. a) The o rbit 1(32) 4 4(2 3) 4 . b) The orb it 12(32) 4 343(23) 4 2. 
c) The orbit 1(32) 4 34(23) 4 2. d) The orbit 1(32) 4 343(23) 4 . 

3. The position on disk 1 of a bounce as a function of the parameter r for some 

periodic orbits in the 4 disk billiard. The bifurcation point is r c = 2.0312 

The positions for r > r c of the two orbits 1(32) 4 4(23 ) 4 and 1(32) 4 343( 23) 4 
are close on the uppermost curve and the two orbits 12(32) 4 343(23) 4 2 and 
12(32) 4 4(23) 4 2 are close on the down-most curve. For r < r c all the (physically 
not admissible) orbits are close. 

4. A sketch of the bifurcation diagram of the Hamiltonian (1). Solid lines cor- 
respond to elliptic orbits and dotted lines to hyperbolic orbits. Only short 
orbits are included. 
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